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Abstract. In recent years various models for the Galactic distribution of massive 
compact halo objects (MACHOs) have been proposed for the interpretation of mi- 
crolensing toward the Large Magellanic Cloud (LMC). A direct way to fit the best 
model is by measuring the lens parameters, which can be obtained by measuring 
the Einstein crossing time and the parallax effect on the microlensing light curve 
and by astrometry of centroids of images. In this work, the theoretical distribu- 
tion of these parameters is obtained for the various power-law Galactic dark halo 
models and MACHO mass functions (MF). For self-lensing as one of the models 
for the interpretation of LMC events, the maximum shift of image centroids and 
the parallax parameter are one order of magnitude smaller than for models with 
dark halos. This can be used as a test for the self-lensing, although the shifting of 
image centroids 0.02 mas, for the LMC events is unlikely to be observed by the 
astrometric missions such as FAME, GAIA and SIM. 
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1. Introduction 

One of the primary candidates for the dark matter in the Galaxy is the massive compact 
halo object (MACHOs) that can be in baryonic or even non-baryonic form. Paczyhski 
(1986) proposed the gravitational microlensing method as an indirect means of detect- 
ing MACHOs. The microlensing effect is the magnification of a background star while a 
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MACHO passes the line of sight of an observer to a background object. Due to the small 
separation between the images produced by the lensing, it is impossible to resolve them 
by ground-based telescopes and only the magnification and its duration can be measured 
by observing the light curve of the background star. 

For one decade many groups contributed in this experiment and observed stars in the 
Large and Small Magellanic Clouds to measure the contribution of MACHOs to the dark 
halo mass and estimate the mass of MACHOs themselves. Due to the degeneracy between 
lens parameters (i.e. mass of the lens, MACHO distance and the transverse velocity), it 
is impossible to obtain them uniquely just by measuring the duration and magnifica- 
tion of events. The degeneracy problem also renders us unable to find the distribution of 
MACHOs in the Galaxy although one can learn much from statistical methods by analyz- 
ing the duration distribution of events. Furthermore, the low number of events observed 
in the direction of LMC causes considerable Poisson fluctuation; to properly constrain 
Galactic models we will need many more events to be detected as will be carried out by 
the next generation of microlensing surveys IjBecker et al. 2004jl . One of the alternative 
ways to overcome the problem of statistical fluctuation is improving the quality of light 
curves, measuring the other observable parameters which can break the degeneracy be- 
tween the lens parameters. 

Apart from the duration of microlensing events, there are two more observable parameters 
that can be measured by accurate observation of light curves. Gould (1992) showed that 
the effect of the rotation of earth around the Sun on the microlensing light curve is a slight 
deviation of the light curve from symmetry. The possibility of observing this (parallax) 
effect has been studied towards the Galactic bulge IjBuchalter fc Kamionkowski 1997jl 
and Large Magellanic Cloud (|Gould 1998(1 . An observational strategy to have sufficiently 
accurate microlensing light curves (with a significant parallax signature) has been pro- 
posed to distinguish two extreme models of the Milky Way and to test the hypothesis 
of self-lensing (Rahvar et al. 2003). However, using the parallax effect is not presently 
sufficient to distinguish between Galactic models with the small number of microlensing 
events. 

The displacement of the image centroids due to microlensing can also be measured 
with high resolution astrometry. The typical angular separations of multiple images 
formed by microlensing in the direction of Magellanic Clouds is 1 milli arcsecond 
(mas) ( |Paczyhski 1996| [Paczyhski 1998| ) and the image centroid is expected to move 
about 0.1 mas. The astrometric missions such as FAME^, SIM^ and GAIA'^ which 
will be launched in 2000s, achieve this high precision. Astrometric accuracies antic- 
ipated for those missions arc: 50 micro arcsecond (/las) for FAME, and 10 iias (or 

^ http://www.usno.navy.mil/FAME/ 

^ http:/ /planetquest.jpl. nasa.gov/SIM/sim_index. htm 

^ http:/ /astro.estec. esa.nl/GAIA/ 
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possibly higher) for SIM and GAIA, where this accuracy depends on the apparent 
magnitude of the stars. SIM can do astrometry for stars as faint as 20th magni- 
tude and it can measure the shift of image centroids of most of the LMC stars 
( |RajagopaI, Boker fc Allen 200 1| ). This positional shift of images can also be used to 
break one more degree of degeneracy. The advantage of astrometric microlensing is that 
the probability of observing a microlensing event is larger than that of photometric mi- 
crolensing ^Miralda,-Escude 19961 m^nik Sahu 2000llHonma 200^. 
In this work we obtain the distributions of event duration, and two parallax and as- 
trometry parameters for various power law dark halo models with various MACHO mass 
functions (MF) and also LMC models (so-called self-lensing) . We show that for the self- 
lensing model, the parallax and the astrometric parameters are one order of magnitude 
smaller than for typical dark halo lensing. 

The organization of the paper is as follows: in Section 2 we discuss the degeneracy prob- 
lem in microlensing and show how astrometry accompanying the parallax measurements 
can break the degeneracy between the lens parameters. In Section 3 we introduce power 
law dark halo models in the Milky Way and obtain the theoretical distributions of event 
duration, parallax and astrometry parameters for various MACHO mass distributions. 
In section 4 we discuss the hypothesis of self-lensing and the distribution of observable 
parameters in this model. The results are discussed in Section 5. 

2. Non-Standard microlensing and degeneracy breaking 

The Einstein crossing time of the microlensing events (te), degenerately depends on the 
mass, distance and the transverse velocity of the lenses. To break this degeneracy, we 
can measure both the deviation of the light curve from the Paczynski form due to the 
parallax effect and the displacement of image centroid during the lensing. In this section 
we give a brief account on these effects. 

2.1. Parallax effect on the microlensing light curve 

As the Earth is rotating around the Sun, its velocity component, projected on the de- 
flector plane, is not negligible compared to the transverse speed of the deflector. The 
apparent velocity of the deflector with respect to the line of sight is a cycloid instead of 
a straight line and this effect causes the light curve of the microlensing event to deviate 
from that of a Heliocentric observer IjAlcock et al. 1995)l . The deviation of the light curve 
from simple microlensing can be described by the parameter 6u which is the projection 
of the Earth orbital radius (a^) in the deflector plane in units of the Einstein radius and 
the projected transverse speed of the deflector {v) on the Earth's orbit. 6u and v are 
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defined as follows: 

5u = a^il- x)/Re, (1) 

where Re is the Einstein radius and x is the ratio of the observer-lens distance to the 
observer-source distance. Equations (1,2) show that by measuring the signature of the 
parallax effect in the microlensing light curves, one can constrain the mass and the 
distance of the lenses. This effect can reduce one degree of degeneracy between the lens 
parameters. 



2.2. Astrometric microlensing 

The angular separations of the images due to the microlensing of the background stars 
at the LMC, lensed by the dark halo MACHOs, are of order of one mas. The image 
centroid sweeps out nearly one m,as during the microlensing event. Not only the ground- 
based telescopes, due to the seeing limit, nor the HST, due to the diffraction limit, can 
achieve the resolution of m,as shifts in image centroid, however observations based on 
interferometry in space may be able to do so. The shift of image centroid during the 
lensing can be derived as follows: 

^ = 2T^^'-' 

where u is the separation between lens and the source in the lens plane normalized to 
the Einstein radius, 6e is the angular size of the Einstein radius {Oe = where D^i 
is the observer-lens distance). For u = \/2 the angular shift obtains its maximum value 
56max = Measuring 56max yields the angular Einstein radius and constrains the mass 
and distance of the lens from the observer. 



2.3. Degeneracy breaking 

The Einstein crossing time, Ie^ depends on three parameters: the mass, M, transverse 
velocity, vt, and distance of lens from the observer, Dg, through 
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On the other hand the parallax and astrometry parameters {6u and 66max) give two 
more equations depending on the mass and the distance of the lens. Using these three 
equations, we can obtain the lens parameters as following: 

(5) 

^sOUOUmax 

1 /c^ap^\ (6) 



a® + \/2DsSuS6r, 
2V2 V G6u 
te Xa^ + y/2DsSu66r, 



a^OUrnax 

""t = -ri — , /Kr. . — • iV 
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3. Constraining Galactic Models via Degeneracy Breaking 
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In the last section, we showed that measuring the duration of events along with the astro- 
metric and parallax parameters can break the degeneracy between the lens parameters. 
In spite of the large optical depth of astrometric microlensing compared to the photo- 
metric one (|Honma 200 l|l , we need to do both astrometric and photometric microlensing 
to get the mass and the distance of the lenses. The mass and distance of lenses can be 
obtained by measuring the two parameters Su and SOmax where their dependence on the 
distance of lenses from the observer is as \J^^- This means that measuring the distance 
of lenses is more accurate the closer the lens is. The fraction of localized lenses among the 
overall microlensing events is expressed by the observation efficiency, which is a function 
of the accuracy of the apparatus and the observational strategy. The expected distribu- 
tions of 5u and 69max can be obtained by convolving the theoretical distributions with 
the observational efficiency of the experiment. 

Here we obtain the theoretical distributions of <£, 5u and 50max in various Galactic mod- 
els for microlensing events in the direction of the LMC. The relevant components of the 
Galaxy, in the LMC microlensing events, are the Galactic disk, halo and the LMC itself. 
For the dark halo, we have considered a large set of axisymmetric models, so-called power- 
law models IjEvans 1994|l . To model the density of the thin and thick disks of the Milky 
Way, we use a double exponential function ( |Binney fc Tremaine 1987| ). Eight power-law 
halo models with a corresponding disk can be combined to build various Galactic models: 
these are the Standard Model (S), Medium halo (A), Large halo (B), Small halo (C), 
Elliptical (D), Maximal disk (E) and two thick disk models (F & G) where we follow the 
designations of Alcock et al. (1996). 

One of the elements for generating i^, 5u and 56 max is the mass function (MF) of 
MACHOs. The tradition has been to use a single MACHO mass (i.e a Dirac-Delta 
function). Here we use two Dirac-Delta and spatially varying MFs. Kerins and Evans 
(1998) and recently Rahvar (2005) have shown that a spatially varying MF for the 
MACHOs, not only can fit the distribution of microlensing event durations in the di- 
rection of LMC, but can also resolve the contradiction between the microlensing results 
and other astrophysical observations such as the unseen white dwarfs in the Galactic halo 
dOppenheimer et al. 200TllTorres et al. 20021 [Spagna et al. 2004| ). The Dirac-Deha MF 
is given by 5{M — Mmacho), where Mmacho in the various Galactic models is given in 
Alcock et al. (1996, 2000). The spatially varying MF has the form MF{r) = 5{M~M{r)) 

I'M \'^/Rhalo 

where M{r) scales as Mu \ 1 , where AIjj and Af^ are the upper and lower lim- 

its of mass of MACHOs, r is the distance from the center of Galaxy and Rhaio is a halo 
scale parameter. The mean mass for the spatially varying MF in the standard dark halo 
model from the likelihood analysis is about the brown dwarf mass (|R,ahvar 2005(1 . The 
physical meaning of this type of spatially varying MF is that the inner part of dark 
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halo contains heavier MACHOs than the outer part. The optimized parameters for the 
spatially varying MF of MACHOs in power-law halo models, to be compatible with the 
LMC microlensing candidates of MACHO experiment is given in Rahvar (2005). The MF 
of the Galactic disk also can be taken from the results of observations with the Hubble 
Space Telescope (Gould, BahcaU & Flynn 1997). 

The distribution of event duration in the various Galactic models can be obtained by 
numerical simulation (Alcock et al. 1996). To obtain the distributions of Su and 69max 
during an observation time of Tabs , we perform a Monte-Carlo simulation, using the prob- 
ability function of lenses along our line of sight. The rate of microlensing events due to 
the MACHOs reside between the distances x and x + dx from the observer is given by: 



where M{x) can be substituted from the desired MF. The location of the selected lenses 
is used to calculate corresponding 6u and 59 max- Figures 1 & 2 show the results for the 
distributions of event duration, parallax parameter and maximum centroid shift in the 
power-law halo models. For the models with dominant halo such as 5, A, B, C and D, the 
distribution of 5u and 59 max looks similar, while they are different from the disk dominant 
models (i?, F and G). The distribution of 59max is sensitive to the type of MFs. For the 
spatially varying MF, the distribution of 59„iax shifts to significantly smaller values than 
single MACHO mass MFs. 

4. Testing the hypothesis of LMC self-lensing 

Due to the low optical depth of microlensing events in the direction of the 
Magellanic Clouds, the EROS and MACHO experiments have obtained only 3 and 
13-17 microlensing candidates, respectively IjLasseree et al. 20001 lAlcock et al. 2000l 
[Tisserand and Milsztajn 2005| ). These studies lead to an upper limit of about 3% (in 
EROS experiment) to 20% (in MACHO experiment) obtained for the contribution of 
MACHOs to the dark halo and a mean MACHO mass of about that of a white dwarf 
IjLasseree et al. 2000l lAlcock et al. 20001 [Tisserand and Milsztajn 2005| ). However, the 
expected white dwarfs have not been seen by studies such as proper motion measure- 
ments | |Oppenheimer et al. 200l]lTorres et al. 20(121 1 Spagna et al. 2004| ). This inconsis- 
tency can be resolved by (i) considering a spatially varying MF for the MACHOs in 
the Galactic halo or (ii) considering the LMC stars playing the role of lenses. If most of 
MACHOs reside in the Galactic halo, the timescale analysis of LMC microlensing can- 
didates implies that the typical mass of lenses is of the order O.SM© while for the SMC 
events, the mass of lenses is estimated to be 2 to 3 Mq ( |Sahu, 2003| ). This inconsistency 
can be solved by considering the lenses as belonging to the LMC. 

The other motivation for the possibility of self-lensing was presented by Sahu (2003), 
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Fig. 1. Distributions of event duration, te, (in days, total number of events is normalized 
to one), parallax parameter, 5u, and astrometry parameter, SOmax, (in mas) for the 
Galactic halo models of S, A, B & C. The solid curves are the distributions derived using 
the Dirac-delta MF and the dashed lines for a spatially varying MF. 



where the number of localized binary lenses was compared with the overall microlensing 
candidates. If half of the MACHOs are considered to be in a binary system (similar to 
the stars), then we expect to observe half of the microlensing events via binary lenses. 
On the other hand only 20% of binary lenses are expected to show the caustic crossing 
feature, then we should see the caustic crossing for 10% of lenses. Two of the Magellanic 
Clouds' microlensing candidates had caustic a crossing feature and by breaking the de- 
generacy between the lens parameters, they are localized within the Magellanic Clouds. 
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Fig. 2. Distributions of event duration, tg, (in days, total number of events is normalized 
to one), parallax parameter, du, and astrometry parameter, SOmax, (in mas) for the 
Galactic halo models of D, E, F & G. The solid curves are the distributions derived using 
the Dirac-delta MF and the dashed lines for a spatially varying MF. 



According to this rough estimation, about 20 microlensing events are expected via LMC 
self-lensing; this number is about the same number of events as have been observed. 
Using the parallax effect and astrometric centroid shift observations can test the hypothe- 
sis of self-lensing. We used a disk model for the LMC ( |Gyuk et ai. 2000| ) and obtained the 
distributions of te-, Su and S0„iax for self-lensing. Comparing Figure Q with Figures 
and (O shows that for self-lensing events du and 59„iax are about one order of magnitude 
smaller than that of halo lensing. For the halo lensing the mean value of the maximum 
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Fig. 3. The distributions of duration of events te (in days), parallax parameter 6u and 
the astrometry parameter 86 max (in mas), for the self-lensing events. The typical values 
of 5u and 69max for the self-lensing events comparing with Figs. (1 & 2) are one order of 
magnitude smaller than for the halo . 



shift for the centroid of images is about 0.2 mas while for the self-lensing, it is about 
0.02 mas. Comparing the resolution of the astrometric missions with the centroid shift 
for the halo and LMC lenses, shows that FAME, SIM and GAIA have enough resolution 
to observe the halo lensing while for the self-lensing events the resolution is of the same 
order as the centroid shift and it is unlikely to the detectable. 
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5. Conclusion 

In this work we study the degeneracy breaking between the lens parameters by measuring 
the duration of events and deviation of light curves due to the parallax effect and the 
maximum displacement of image centroids during the lensing. 

We used various Galactic models and obtained the theoretical distributions of observable 
parameters. In order to have the observed distribution of parameters we need to convolve 
these distributions with details of the observing program such as the photometric and 
astrometric accuracy and the sampling rate. One of the applications of astromctry with 
the parallax measurements is testing the hypothesis of self-lensing. The maximum shift 
of the centroid of images and parallax parameter for the halo lensing are one order of 
magnitude large than the self-lensing ones. 

In practice, observations to measure the three parameters of microlensing can be done 
by a ground-based telescope accompanying a space-based interferometry telescope. 
The ground-based telescope alerts one to the microlensing events and undertakes 
photometric measurements and the space-based telescope measures the displacement of 
image centroids. Rahvar et al.(2003) proposed a strategy for observation of microlensing 
events based on both alert and follow-up telescopes. They obtained the observational 
efficiency for the detection of parallax effect and estimated the number of microlensing 
events essential to distinguish between two extreme Galactic models and to test the 
self-lensing hypothesis. This type of observational strategy can be extended by adding a 
space-based astrometric telescope. 
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